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R ecently I was working through some problem sets on determining 
volumes by triple integrals in cylindrical and spherical coordinate 
systems. The textbook I was using included many interesting problems involv- 
ing spheres, cylinders and cones and the increasingly complex solids that 
arose from the intersections of two or more of these. I was surprised it did not 
include the common torus. So naturally I tried it to see what does happen. 

I am still unsure why the textbook and several others I have consulted do 
not include the torus in the problem sets for triple integration in cylindrical 
coordinates. The technique produces the solution quite neatly. However, 
after many hours and pages of working, I quite understand why problem sets 
on triple integrals in spherical coordinates avoid the torus. It is a long and 
arduous journey. It is however a very interesting journey with many unex- 
pected sights along the way. I have written this brief guide for readers who 
need some assistance with this trek. Pack some doughnuts for visual aids and 
sustenance! 

Of course, there are much simpler ways to determine the volume of a 
torus, which raises the question: why are mathematicians willing to spend 
hours solving a problem by a difficult method when they already know an easy 
method? 

Rather than wasting time in lengthy argument on that issue, I will simply 
offer my three favourite answers: 

1. because I just do not like to be beaten by a problem; 

2. because that behaviour is inherent in the definition of “mathemati- 
cian”; 

3. because it is fun. 

Introducing the torus 

Consider a circle in the xy-plane with centre {R,0) and radius a< R This is the 
circle 

{x - R)^ + J 
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Rotate the circle around the 31-axis. The resulting solid of revolution is a 
torus. It is sometimes described as the torus with inner radius R - a and outer 
radius i? -1- a. It is more common to use the pro numeral r instead of a, but later 
I will be using cylindrical coordinates, so I will need to save the symbol rfor 
use there. 

Before considering cylindrical and spherical coordinates, it is useful to 
briefly review three simple well-known methods of determining the volume of 
this torus. 


Pappus’ Centroid Theorem 

The shortest method is to employ Pappus ’ Centroid Theorem. 

The circle of radius a has area Tzc^. When it is rotated around the 31-axis its 
centroid, which is simply its centre, tracks a circle of radius R. Thus the centre 
travels a distance of ‘J.nK Pappus’ Centroid Theorem then immediately gives 
the volume generated as 


V= na^ X 2 nR = 2 n‘^a^R 

For this method to be convincing to students they need to prove Pappus’ 
Centroid Theorem, but the proof is within the reach of students who are 
studying triple integrals in cylindrical or spherical coordinates. 


Two slicing techniques 

There are two “volume by slicing” techniques that allow the result to be 
readily determined with a single integral. These will only be described very 
briefly since they are both well known. For example, in New South Wales, 
Mathematics Extension 2 students should encounter both these methods in 
Year 12 . 

In both cases, rather than rotating the full circle {x - R)'^ + around 

the 31-axis to produce the torus, it is simpler to use only the semicircle above 
the x-axis to generate the “top” half of the torus and double the result. 

Take a thin horizontal slice of the semicircle between 3 and 3 -(- 83. Rotating 
this around the 3-axis gives a prism of an annulus. I t has h eight 83. The 
ann ulus at t he base of this prism ha s inner radius R — \Ja^ - 3^ , outer radius 
R+^a^ - 3^ and thus area AnRyja"^ -3^ . The 
required volume is 

a a 

y = 2 j 4 :nR-^a^ -/^dy = 8 nR j -^a^-y^dy 
0 0 

While the final integral can be evaluated by using a trigonometrical substi- 
tution, it is simpler to recognise it as the area of a quadrant of radius a, giving 

V = 8 nRx — na^ = 2 n^a^R 
4 
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The other “volume by slicing” method involves taking a thin vertical slice 
of the semicircle between x and x + 5x. Rotating this around the y-axis 
produces what is often called a cylindric al sheik thou gh it is reallyjust another 
prism of an annulus. It has height y = —{x — Rf . The annulus at the base 

has inner radius x, outer radius x + 5x and thus area 27ix 5x, ignoring the 
second order terms in 5x. The required volume is 

R^d I 

V = 2 j 2nx'<Ja^ -(x-R)^ dx 

R-a 

R-\-a . 

= 47t [ xyja^ —{x — R)^ dx 

R-a 


The substitution u = x - R produces 

a j 

y = 471 J {M + R)^ja^ —u^du 


a p 

J Wc 


= 471 I u'^ du + AnR I '^a‘^—u‘^du 


a ^ 

J ^^C 


2 .,2 . 


= 471 


1 / , 
3 




+ 47ti?x — Tia^ 
2 


= 2n^a‘^R 


where the second integral has been evaluated by recognising it as the area of 
a semicircle of radius a. 


Reorienting the torus 

Cylindrical and spherical coordinate systems often allow very neat solutions to 
volume problems if the solid has continuous rotational symmetry around the 
z-axis. While I will use the same torus as discussed above, it will be oriented 
differently relative to the axes to produce the required symmetry. 

Readers who can be trusted with sharp knives should now take a conven- 
ient doughnut and place it on a suitable clean horizontal surface. Use a sharp 
knife to slice it in half with a plane parallel to the horizontal surface. Slide the 
bottom half out and dispose of it in the obvious manner. (You will need to 
keep your strength up to cope with the algebra that follows!) The newly cut 
annular surface of the doughnut should now be sitting on the horizontal 
surface. Place the origin on the horizontal surface at the centre of the hole in 
the middle of the doughnut. The horizontal surface is the xy-plane with the z- 
axis rising vertically from that plane. I will determine the volume of the 
half-doughnut and double the result. 

Readers who do not have a doughnut to hand will have to content them- 
selves with taking the circle {x-R}^ + z' = c^, y= 0 and rotating it around the 
z-axis to produce the required torus. Consider the half of the torus above the 
xy-plane. 
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Cylindrical coordinates 


In cylindrical coordinates, a point is located by the triple (r, 6, z) where z is 
the usual rectangular z-coordinate and (r,9) are polar coordinates in the xy- 
plane, 6 being measured anticlockwise from the positive x-axis. For an 
arbitrary 9 draw an r-axis in the xy-plane at an angle 9 anticlockwise from the 
positive x-axis. The cross section of the torus in the rz-plane is as follows 
(Figure 1). 


z 



Figure 1 


This cross-section is the same for all values of 9. This property hints that a 
solution by cylindrical coordinates is likely to be efficient. Incidentally, it also 
means the above figure looks identical to the cross section in the xz-plane. 

When drawing the cross-section, only consider positive values of r. 
Allowing 9 to run from 9 to 271 generates the whole torus, so there is no need 
to consider negative values of r. 

The circle shown above is {r- R)^ + ^ = a^. The top half of the circle can 
be generate d by allowin g r to run from R — a to R + a and allowing z to run 
from 9 to '\ja^ —{r — R )^ . Thus, using a triple integral in cylindrical coordi- 
nates the volume of the torus is 

271 R+a -^a^-ir-Rf 

y = 2 [ [ { dzr dr dQ 

e = 0 r=R-a z=0 

It was noted above that the cross section was independent of 9. As a result 
of this the inner two integrals are constant with respect to 9, and so they can 
be taken outside the outer integral as a common factor, giving 

271 R+a -^a^-ir-RF 

V = 2 \dQ f \dzr dr 

e = 0 r=R-a z=0 

R+a 

= 2x27t J rxj - (r - R)^ dr 

r=R-a 
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This integral, with the dummy variable r replaced by x, has already been 
evaluated in the last of the simpler methods given above, the result again 
being 

y= 2n^a^R 


Spherical coordinates 

In spherical coordinates a point is described by the triple (p, 0, (|)) where p is 
the distance from the origin, (|) is the angle of declination from the positive z- 
axis and 9 is the second polar coordinate of the projection of the point onto 
the xy-plane. 

Allow 6 to run from 0 to 271. For any 9, the cross-section in the resulting rz- 
plane is as shown in Figure 2. Only consider the half of the torus above the 
xy-plane. Let OB be a tangent to the semicircle. 

z 



For the arbitrary 9, determine the integration limits for (|). Imagine (|) as 
controlling a ray in the rz-plane, the ray being able to pivot around the origin. 
Start with (|) = 9, meaning the ray overlaps with the positive z-axis. Allow (|) to 
increase, so the ray rotates clockwise around the origin. The ray first intersects 
the semicircle when it overlaps the tangent OB, where (j) = (3. Triangle OAB is 
right-angled, so 

cos P = — 

R 


The ray continues to intersect the semicircle until the ray overlaps the r-axis 
where 



Hence allow (|) to run from p 


= cos 


r \ 
a 


to 


n 

2 ' 


Now consider an arbitrary (|) within this range and determine the integration 
limits for p (Figure 3). 
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z 



The points to be included in the integration are those making up the inter- 
val CD. That is, the minimum acceptable value of p is the length of OC and 
the maximum is the length of OD. The lengths of OC and OD can be deter- 
mined by applying the cosine rule to triangles OAC and OAD respectively. Let 
j denote the length of OC in the first case and OD in the second. Curiously, 
both triangles give the same equation. 


a^ = f+R^ 
= f+R^ 


-2 jR cos 


f 

n 


12 


y 


-2ji?sin(|) 


This is a quadratic equation in j. 

f -{2Rsm(^)j + R^ -a^ =0 


Solving gives 

2R siml) +J4i?^sin>-4(i?^-a^) 

7 = 

2 

= R sincj) sin^(t)-|/?^ -a^j 


or 

J = R sin(]) —R^ cos^cj) 

That last rearrangement is probably not an obvious thing to do and it is 
not clear whether it improves things or not. If the rearrangement is not done 
now, then at a later point in the working it will become the obvious thing to 
do. I will do it now since it does marginally simplify the next few steps. If this 
were being set as a problem for students, th is is a good po int to provide guid- 
ance by asking them to “Show j = R sin(|) inja ^ -R^ cos^tj) .” 

The two solutions to this quadratic are the lower and upper integration 
limits for p, corresponding to the cases where J is the length of OC and OD 
respectively. It now becomes important to distinguish between the two solu- 
tions. Let the lengths of OC and OD be ji and 72 respectively. That is. 
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= J?sin(|) —R^ cos^(|) 

= R sin(j) 

After lengthy algebra such as the above some reasonableness checks are in 
order. 

Consider the case (j) = (3. This makes the points C and D coincide at the 
point B of Figure 2. That is, the minimum and ma ximum acceptable values of 
p should coincide at the length of OB, which is —a^ . Substituting (|) = (3 
into the expressions for and gives 

= R sin (3 — \ja‘^ -R^ cos^ P 

^2 = R sin (3 +^ja^ -R^ cos^ (3 


where cos B = — . Hence 
R 


—R^ cos^ (3 = 


a^-R^ 


( ^2 
a 


= 0 


71 


Also, since 0 < B < — , 

2 


Hence 


sin [3 = -^1-cos^ (3 = 


sIr 




R 


h=R 


-^R^-c 
R 

=4r^ 

4r^c 


J2 ~ R 


R 


= -\Ir 




as expected. 


The other easy check is (j) = — . 

In Figure 3 this places points C and D at the intersections of the semicircle 

with the r-axis. The resulting lengths of OC and OD should be i? ± a. 
7Z 

Substituting ^ the expressions for ji and gives 

j, = Rsin — --, -R^ cos^ — 

2 V 2 

= R — a 

= Rsin^+^a^ -R^ cos^ ^ 

= R + o 

as required; that is, both checks produced reasonable results. 

Recall that j] and j 2 are the integration limits for p. Thus, using a triple 
integral in spherical coordinates the volume of the torus is 
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271 2 J‘i 

V = 2^ j Jp^sin(|)^^pd(l)de 

e=o <|)=Pp=ji 

Since the cross section is independent of 0 the inner two integrals are 
constant with respect to 9, and so they can be taken outside the outer integral 
as a common factor, giving 

71 

2n 2 J 2 

V = 2 \dB \ [ p^sin(|)dpd(|) 

e=0 P = Pp = ;i 

K 

= 2 x 271 j ^(jl-jfjsin^# 

71 

4 2 

" q" 1 (i2-h)sinc^# 

p=p 

The next step is to evaluate jl - jf where 

= R sin(|) —R^ cos^(|) 

^2 = R sin(|) +-^( 2 ^ —R^ cos^(|) 

and to hope that the product of this result and sin (|) is integrable. 

This step is ugly howev er it is attemp ted. Perhaps the least ugly method is 
to set u = Rsin ^ and v= -yja^ -R^ cos^ ^ , so thatji = u-v andj 2 = u + v. Then 
wade into the algebra and find 

il - if = (i2 - ii){i2 + hh + if} 

= +2uv + v‘^^ + [u^ —v^^ + {u^ —2uv + v‘^^ 

= 2u{3M^+n^| 

= 2 -^ 0 ^ —R^ cos^(|)|3i?^ sin^(|) + a^ —R^ cos^{|)| 

= 2Va^-f?^cos>{a^ + 3/?^ - cos^ (|)| 

By this point any sane person would have concluded that spherical coor- 
dinates are just not the bright way to solve this problem, but people who use 
mathematics for recreational purposes are made of sterner stuff and will press 
on. 

Substituting this result into the expression for the volume gives 

TT 

y = — 7t j -R^ cos^(|)|a^ -l-3i?^ -4i?^ cos^(|)|sin(|)d(t) 

This integral does not look very promising. However, the volume of a torus 
does have a nice solution so further effort is warranted. The obvious line of 
attack is to use a substitution to tidy up the mess within the square root sign. 
Hopefully some unexpected side effect of this will deal with the troublesome 
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sin (|) factor and allow the integral to be evaluated before I run out of differ- 
ent letters to use as pronumerals. 

Try the substitution R cos (|) = a cos w. 

Hence cos^ (|) = «^ cos^ w and R sin di^ = a sin w dw. Turning to the 
limits: 


([) = P = cos 




yRj 


w = 0 and (b = — ^ w = — 
2 2 


The substitution has simplified the lower limit, an unexpected bonus. 
Applying the substitution gives 

K 

^ P I 2 2 2 r 2 2 2 2 1 ^ 

y = — 71 J -a cos w^a +3i? -4a cos rcj — %\nw dw 


n 

2 


= — n I" asinwja^ -4a^ cos^ wj-^sinw dw 

J I iR 


w=0 


8a TZ 
8R 


2„ 2, 


r 9r9 9 9 91 

sin wia +SR -4a cos w>dw 


Finally, the integral is looking manageable. Using the results 

sin^ w = i (l - cos 2w) 

and cos^ w = ^(l-l-cos2w) 


gives 


8a^n f 1 


V = J — (l-cos2w)|a^ -t 3/?^ -2a^(l-tcos2w)|dw 

w^O ^ 

TC 

2 

SR 2w){3i?^-a^-2a^ cos2wjdw 


SR J 2 

td=0 

TC 

4a^7t ^ 


JT=0 


4a^7i 


SR 


j SR'^ -a^ -(SR^ +a^)cos2w-t2a^ cos^ 2 

v=0 


w dw 


Then cos^ 2w = i^l + cos4w) gives 


T = 


4a^7t 


3i? 


1 SR^ —a^ —^SR^ jcos2w-l-a^(l-l-cos4w)dw 


zt;=0 

n 


A ^ ^ 

^ J 3i?^ — ^3i?^ +a^ jcos2w + a^ cos4w dw 
3-/? ^ 

i!t;=0 


Symmetry could now be used to argue that the terms in cos 2w and cos do* will 


57 


The volume of a torus using oylindrioal and spherioal ooordinates Australian Senior Mathematios Journal 19 [2] 



Australian Senior Mathematics Journal 19 (2] Farmer 


evaluate to zero on integration over 0 to g . Alternatively, integrating in the 
normal manner gives 


V = 


4a^K 


3R 

4a^n 




1 


?>R w 3i? -\-a \sm2w-\- — a sin4rt» 


If 


-<0 


3R 
= 2n'^a'^R 


X 3R^- 


Of course, I am not suggesting that using triple integrals in spherical coor- 
dinates is an efficient way to find the volume of a torus. Given the great 
potential for algebraic errors, it does not even make a particularly good exam 
question. It is more the type of problem that could be tackled as a group exer- 
cise, where students can check their progress with each other frequently and 
find the errors before they carry them through too many steps. 

What I find most fascinating about this problem is the number of different 
mathematical concepts that need to be applied to reach the result. Students 
sometimes seem to treat mathematics as a series of unrelated techniques each 
with their own clearly demarcated area of application. This problem demon- 
strates the need to reach for tools that would not normally be thought of as 
part of the area of knowledge labelled “triple integrals”. 
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